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In calculation of the theoretical scattering coefficient for a terrain, previous authors 
tentatively assumed the normalized autocovariance function p{r)=e~'^^' for the ground 
elevation as a function of distance from a given point. Recently autocorrelation studies 
were made using maps with contours ranging from one to twenty-five feet. These n^snlted 
in curves of p{r), w^hich are approximated by exp {—\r\/B). The theoretical scattering 
cross section (o-q) of many such terrains can be expressed as 
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where a-, X, k, and $ are standard d(^viation of the target terrain, wavelength, wavenumber 
(27r/X) and the angle of incidence respectively. For the case where \/B is small as compared 
to k, the above expression becomes 
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These expressions, wlien normalized, are in agreement witli experimental results of oti'er 
authors. It is also noteworthy that the results obtained with an acoustic sinudator model 
compared very well with this theorc^tical expression. This work is based on the i:roperty 
that the ground is conducting and has random elevation variations. Theoretical results 
calculated on the basis of varying ground impedance rather than its elevation are also in 
agreement with this expression. 



1. Introduction 

In recent years, the calculation of backseat ter 
from a rough surface, with very obvious exterision 
and application to tlie radar return from the moon, 
has attracted considerable attention. Here we 
follow the basic approacli of Davies [1954] as modi- 
fied by Moore [1957] and (^ooper [1958]. The 
modified Kirchhoff-Huygens' princi])le is employed 
in the calculation using modified s]))ierical variables 
of integration. The ground model used assumes 
'^ facets,^' of variable size and height above mean 
ground level, whose position is described statistically. 

Instead of assuming a correlation model for 
terrain, as is usually done, tlie present approach 
used contour maps of different terrain samples in 
the United States to calculate tbe terrain-elevation 
autocovariance and other statistical ])i-operties. 
The overall approacli to the problem is an approxi- 
mation, but t Jie results so obtained are very reason- 
able indeed, insofar as the comparison with terrain 
return and moon-echo data reported in various 
pubhcations [Briggs, 1960; Hughes, 1960; Nielson, 
1960] is concerned. 



1 This work was sponsored bv tho Naval Ordnance Test Station, China Lake, 
Calif., under contract No. N 123 (60530) 181 38 A. 



2. Statistical Properties of Terrain 

Various autocovariance [Anient, 1953; Davies, 
1954; Moore, 1957b; Daniels, 1960] functions for 
rough terrain have so far been assumed. In some 
instances [(^ohen, 1948] an effort was made to 
derive an approximate expression for the autocorrela- 
tion from actual (hit a on terrain elevation. A 
general expression [Norton, 1960] for space correla- 
tion function p (?) describing the random variation 
of the refractive index over space is given here: 



where 
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/q= characteristic scale 
^1= constant 
r(/x)=gamma function 
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modified Bessel function of the second kind. 



It seems to cover some of the most commonly as- 
sumed [Wheelon, 1959] expressions for correlation 
functions. 

A search of the literature shows that small-scale 
perturbations of the terrain elevation have not so far 
been used to calculate an experimental autocovari- 



427 



Table 1. 



Location of terrain samples and the surface 
characteristic constant, B 



Sample 


General location 


Average 
elevation 


Standard 
deviation 


Sampling 
interval 


Surface 
characteristic 
constant, (1/B) 


1 

2 ._- 

3 - 


Lawrence, Kans 

Tees Nos, Ariz 

Turtle Mountain, 
N. Dak 


n 

853. 5 
5254 

1990 
2617 

5028 

7530 
6780 


ft 
5.63 
49 

8.2 
13.6 

5. 18 

1030 
250 


ft 
30 
5 

5 
5 
5 

62.5 

20 


0. 1193 X 10-3 
1.3213x10-3 

0. 7070 x 10-3 


4 

5 

6 


Gila River, Ariz 

White River, Ariz__ 
Mountain Park, 
N Mex 


. 9236 X 10-3 
1.3610X 10-3 

0. 2729 X 10-3 


7 


Sandia Park. 
N. Mex._ .___ 


. 9860 X 10-3 









ance function. For this reason, contour maps for 
seven different terrains in the United States were 
selected, as listed in table 1. This set included 
relatively flat land, rolling plains and some rugged 
mountain areas. The average elevation (above sea 
level) and standard deviation of these samples varied 
from 850 to 7,530 ft and 5 to 1,030 ft, respectively. 
Three distinctly different but random lines were 
drawn on each of these contour maps. Along these 
lines, the elevation of the terrain was read to within 
one tenth of a foot (by interpolation) on 1-ft contour 
maps, and within 5 ft on 20- and 25-ft contour maps. 
Using a reasonable compromise between resolution 
and confidence level [Blackman and Tukey, 1958], 
and due to the unavailability of less-than-one-foot 
contour maps, horizontal sampling intervals of 5, 20, 
30, and 62.5 ft were used. The number of points 
for a subsample varied from 58 to 610. The auto- 
covariance was calculated using a CRC-102 com- 
puter. The resulting curves were averaged for each 
type of terrain. A simple theoretical approximation 
of these curves for high confidence level portions was 
found to be 



where 



.{r): 



o-\r\IB 



(2.2) 



B 



characteristic constant 
r= distance between points 



Seven cases for the experimental autocovariance are 
shown in figure 1. 

In figure 1 the autocorrelation curves for samples 

1 and 3 seem to slope off faster than the rest. In case 
of sample 1 this is caused by the lag distance for a 
given number of lags being greater tban for samples 

2 through 5 and 7. In view of this, if all these 
curves were plotted on a semilog paper (as in fig. 1), 
for a given lag distance rather than for a fixed number 
of lags, these would plot approximately as straight 
lines over limited distances. These two curves can 
also be approximated as p(v) = e~^'^^'^ cos cr, where c 
was found to be 0.2188 and 0.7931 deg/ft for samples 
1 and 3, respectively. For near vertical incidence, 
the distances involved are rather short, and therefore 
cos cr is approximately unity. It is felt that the 
apparent cosine factor in the autocovariance appears 
because of the relatively short length of the sample 
used in this case. This seems to indicate that the 
autocovariance function p(r) often varies as e"''"'''^. 
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Figure 1. Normalized autocovariance curves — {portions with 
highest confidence level). 



Sample No. and type of terrain 



1 (relatively flat) 

2 (rolling and flat) _ _ . 

3 (rolling). _. 

4 (irregular, sloping). 

5 (flat, slight slope) _ . 

6 (very rugged mtn.) 

7 (rugged mountain) . 



Length of data 

(points in each 

subsample) 



58, 101, 079 
252, 276, 215 
189, 194, 144 
226, 196, 261 
287, 264, 272 
216, 187, 187 
187, 187, 187 



Sampling 
interval 



30 

5 

5 

5 

5 
62.5 
20 



It is known [Cohen, 1948; Isakovitch, 1952] that 
the slope of p(r) at r equal to zero should be zero 
unless the ground elevation function has infinite 
slopes. In this instance, p{r) has a slope of —1/B 
which is in most cases believed to be of the order of 
10~^ (see table 1). Considering the minute size (in 
the general range of 10"^ to 10~* m) of the earth parti- 
cles, this may be quite a reasonable approximation to 
the exact description of the terrain roughness. 

It was also found that the ground elevations 
generally are normally distributed. For a certain 
type of terrain, the ground elevation data along a 
random line is most probably a random function 
belonging to a large, approximately stationary en- 
semble. It is further suspected that the expressions 
obtained for the scales of roughness covered in this 
study when reduced by a scale factor, might well 
describe extremely minute variations of ground 
roughness as seen by very high frequency waves. 
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3. Scattering Coefficient 

The expressions for vector waves, in radar return, 
althoiit^li easy to set up are very difficult to evaluate. 
For this reason the Kirchhoff-Huygens' principle has 
been applied to scalar waves. This leads to difficult 
inl(^gration problems and other such complexities, so 
certain simplifying assumptions were made. Similar 
assumptions were made by others [Davies, 1954; 
Moore, 1957b]. These are as follows: 

(1) No portion of the ground is shielded from the 
incident radiation. 

(2) The ground is considered to be a perfect 
conductor. 

(3) The magnitudes of the surface currents are of 
the same order as those of a plane reflector, but the 
phase varies in a random manner, depending on the 
height of a particular point. 

(4) The reradiation from a particular small area 
on the ground is isotropic. 

(5) The antenna gain G is essentially uniform for 
— ^o<C^<C^o, and is zero outside of this range. 

(6) The results ol" tliis scalar wave approach are 
an approximation to those for the vector waves, as 
most field vectors would be nearly parallel to the 
surface involved 

The variables </> and 0' were replaced with yp and 
^', where yp^ is the sum of yj/ and a as shown in figure 2. 
This change of variable was made following Davies' 
[1954] work, in order to make the integration a little 
less complex. 
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Figure 2. 



co'o=a' 



AO'B=v|/ 
AO'C=v|/' 



Geometry of the 'problem ^ showing interrelation 
between various variables of integration. 



The radar equation as applied to pulse radar 
[Moore and Williams, 1957a], gives the average re- 
ceived power Pr from many scatterers as 
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Pr = 



where 
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dR. (3.1) 



Py— power transmitted 
7^= range 
0-0= scattering coefficient per unit area (assumed 

independent of 0) 
^=:time 

c= velocity of propagation (velocity of light) 
G^=: antenna gain (assumed independent of <t>) 
X= wavelength 
d, </) = angles as shown on figure 2. 

The field E at 0' is obtained by applying Huygens' 
principle to the nullified geometry of figure 2, and 
can be written as 

E= r^ (7,77)1/2 ^os de-^'^'^'HA. (3.2) 

where 

77= intrinsic unpedance of free space 

/ d\p \ 

^M=area element on the grouml l=RdR ^j^ )* 

One arrives at the following expression for E by 
applying (3.2) to a rough surface. 



CT 

If^ f^o cot e / PrGy y^ 

xj-cj-^o R Wr'J 



exp [j2k(R-d cos d)WdR (3.3) 
where 

Jc=wsive number ( -^) 

r= pulse width 
2^0= beam width of the antenna 
5=6(jBo+s,iA+o:)=elevation of ground height 
above mean ground level at a point located 
at range Rq+s and modified azimuth angle 

Bo=range to the projection of a point at range R, 
on the mean ground level (see fig. 3). 

Equation (3.3) takes into account the variations in 
phase and neglects the other effects of the changes in 
range. This is based on the approximation that the 
percentage changes in range are not appreciable for 
near-vertical incidence. The received power is 



F,=^„M^'^^ 
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(3.4) 
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MEAN GROUND LEVEL 

Figure 3. A typical point, at a random height 5o above mean 
ground level. 

where E^ stands for complex conjugate of E and i^e 
is the symbol for ^real part of ^ the expression follow- 
ing it. The difference in ranges of two points on 
ground located at (Rfi^yp) and {R^,d\\[/^) can be ex- 
pressed as follows with the help of figure 3, 

B-R' = [Ro-8(Ro,rp) cos e]-[R'o-8(Ro,rp') cos 6'] 

= —s+[K^o+s,\p+a) cose' — 8(Ro,\P) cos 9] 

= -s+(d—do) cos (9 (3.5) 

where ^?=^^' (assumed) 

8=8(Ro+s,xl^+a) 

8o=8(Ro,rP) 

Ro^Ro+S 

The substitution of (3.3) in (3.4) results in 

1 ^ CPtG^ cos e cos 6' 



^'~32w' ^^ 



P 



R^R^' 



(3.6) 



where the primed quantities refer to a point at (R^, 

e', ^'). 

The average received power can now be expressed 
in terms of R, Sj xp, and a by substituting (3.5) in 
(3.6) as 



327rV _^J _20^J _ (RQ-h)JEQ 



^0 Pj.G^ cot e cot 6' 



RR' 



exp [j2k{s+(8o-8) cos d]]dRdsd^da (3.7) 

where primed quantities refer to values at a point 
different from that for unprimed quantities. The 
limits —2^0 to 26^ on the variable yp are used to cover 
the entire 360° of the azimuth angle. 



It is apparent from the statistical studies of contour 
maps that the probability density function for ter- 
rain elevation above the mean (say 8) can be written 
as 



p{5)- 



1 



(2ir<r2)' 



-«2/(2<r2). 



(3.8) 



The conditional (normal) density function for the 
elevations 5 given 8q at two points on the ground, a 
distance r apart, is given by 



p{8-M,r) = 



[27r(72(l-p2)]i/2 



e 2<T2(l-p2). 



(3.9) 



The expression exp [j2k[s+{8Q—8) cos 6]] in (3.7) 
is first averaged over 8 using (3.8) and the result is 



exp [j2ks—a^¥ cos^ d+j2k8o(l — p) cos 6], 



(3.10) 



where 



--2a'(l-p'), 



Now (3.10) is averaged over 8o using (3.9), and it 
becomes, 



exp [j2ks—^kV{l-p) cos^ 6], 



(3.11 



One can deduce from Davies' [1954] approximation 
r^^RW+s^ cosec^ 6 (see fig. 2) that 



r^Ra- 



s^ cosec^ d 
2R^ 



(3.12) 



It can be now shown that the integration of (3.7) 
after substituting (3.11), (3.12), and (2.2) results in 
an expression, which when rearranged, becomes 
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j\ R3 J B\^ Vsin d) 
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BU J 



dR. (3.13) 



Here it has again been assumed that R^ and 6^ are 
approximately equal to R and 6 respectively, for 
near-vertical incidence. From the comparison of 
(3.1) with (3.13) it is clear that the scattering coef- 
ficient, (To, is given by 



(To= 



4V27r52 



Vsm 6/ 



p — ik^a^ cos 2 d 



1 (ri.-l)l [2P52 sin" e+n?f'^' ^"^^ ^ 
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Tlio value of \/B is far less tluiu k for nearly smooth 
surfaces, and. lieuce the scattering coefficient for such 
surfaces can be approximated from (3.14) as 

4o-2 
(To (for nearly smooth surfaces) — -ttt (0 cot^^) (3.15) 

4. Experimental Verification 

For nearly smooth surfaces, the surface character- 
istic constant 1/B^O, and (3.15) gives the scattering 
coefficient. This result compares very closely with 
published results [Nielsen, 1960] for new ice as shown 
in table 2. 

For rough (not nearly smooth) surfaces, (3.14) 
describes the relationship of tlie scattering coefficient 
(To and other variables such as the angle of incidence 
6, wavelength X, standard deviation a and surface 
CO variance constant B, etc. Two curves of the scat- 
tering coefficient ctq versus 6 for each of the three 
values of X/7^, 0.1, 0.5, and 1.0 for a/\ equal to 0.05, 
and 0.1 are shown in figure 4. Tt may be noticed 
that as the surface beconu^s rougher, or as \/B 
increases for a specified X, the scattering coefficient 
curve becomes flatter, showing t\w relative impor- 
tance of the contribution of the power return from 
the surface at angles other than tliose near zero. 
As expected, when the surface becomes smoother or 
1/B decreases, the received power seems to come 
primarily from near-zero angles. These curves are 
quite similar to those recentlv published [Campbell, 
1959; Dye, 1959; E(Hson, 19()()]. Tlie experimental 
data [Nielson, 1960] on desert and new ice also seems 
to follow the pattern of these theoretical curves de- 
scribed above. 



Table 2. 



Comparison of theoretical versus experimental scat- 
tering coefficient 
(jSTormalized) 



00 


(TO Theoretical 


<ro 


Experimental 


30 


1.000 




1.000 


40 


0. 291 




0.308 


50 


.088 




.089 


60 


.022 




.021 


70 


.004 




.004 



Similar results ^ obtained by an acoustic simulator 
at the University of New Mexico also verify tliese 
theoretical conclusions. It is interesting to note 
that an experhnental expression {(To^^ (Tie~^^^ for 
6 in radians) for the scattering coefficient [Hughes, 
1960] of the Moon, verified by photographic astro- 
nomical [Briggs, 1960] calculations, is a very close 
fit to the graph of eq (3.14) for o-/X=0.1 and X/5=1.0 
in the range of incidence angles of 3 to 14 deg. The 
authors believe that other values of B and cr may 
well be appropriate, as only a limited set has been 
tried to date. 



2 A. R. Edison and R. K. Moore, Preliminary report on an acoustic simulator 
for investigation of backscatter of E.M. waves, unpublished report of the Univer- 
sity of New Mexico, Albuquerque, N. Mex. (1961). 
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Figure 4. Scattering coefficient versus the angle of incidence. 



5. Conclusions 

The scattering coefficient (o-q) for nearly smooth 
surfaces is inversely proportional to the wavelength, 
but varies- directly witli (a^), {6 cot^ 6) and l/B^ where 
(T, 6, B are standard deviation, angle of incidence, 
and the terrain characteristic constant respectively. 
For rough surfaces it has a negative exponential 

factor, where the exponent is made up of r-^ — 

times a constant. The surface characteristic con- 
stants B and (t can be calculated from the radar 
return data. Although approximate, the theo- 
retical results agree well with the experimental data; 
and therefore, suggest the usefulness of the approach. 
Tlie application of these results may be extended to 
the moon-echo data, with proper corrections for 
Faraday and liberation effects, etc. This investi- 
gation has established that for near-vertical inci- 
dence, the normalized autocovariance for the terrain 
elevation is more often of the exponential form exp 
(— \r\/B) rather than the Gaussian form, exp (—r^/B), 
The former may well be more appropriate for finer 
terrain irregularities than those considered in this 
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paper. It may also be representative of the normal- 
ized autocovariancc of the moon surface. An exact 
theoretical, but usable, expression for the scattering 
coefficient has been very reasonably approximated 
by the results of this study. 
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